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While slowly turning the ends of a single molecule of DNA at constant applied force, a discontinuity 
was recently observed at the supercoiling transition, when a small plectoneme is suddenly formed. 
This can be understood as an abrupt transition into a state in which stretched and plectonemic DNA 
coexist. Wc argue that there should be discontinuities in both the extension and the torque at the 
transition, and provide experimental evidence for both. To predict the sizes of these discontimiitics 
and how they change with the overall length of DNA, wc organize a phcnomenological theory for the 
coexisting plectonemic state in terms of four parameters. We also test supercoiling theories, including 
our own elastic rod simulation, finding discrepancies with experiment that can be understood in 
terms of the four coexisting state parameters. 



A DNA molecule, when overtwisted, can form a plec- 
toneme [1, 2] (inset of Fig. 1), a twisted supercoil struc- 
ture familiar from phone cords and water hoses, which 
stores added turns (linking number) as 'writhe.' The 
plectoneme is not formed when the twisted DNA goes 
unstable (as in water hoses [3]), but in equilibrium when 
the free energies cross — this was vividly illustrated by 
a recent experiment [4] (Fig. 2), which showed repeated 
transitions between the straight "stretched state" (SS, 
described by the worm-like chain model [5] ) , and a coex- 
isting state (CS) of stretched DNA and plectoneme [6]. 
This transition, in addition to being both appealing and 
biologically important, provides an unusual opportunity 
for testing continuum theories of coexisting states. Can 
wc use the well-established continuum theories of DNA 
elasticity to explain the newly discovered [4] jumps in 
behavior at the transition? 

The recent experiment measures the extension (end- 
to-end distance) and torque of a single molecule of DNA 
held at constant force as it is slowly twisted [4]. A 
straightforward numerical implementation of the clastic 
rod model [7-9] for DNA in these conditions (with fluc- 
tuations incorporated via entropic repulsion [7]) leads to 
two quantitative predictions that arc at variance with 
the experiment. First, the experiment showed a jump 
Az in the extension as the plectoneme formed (Fig. 1) 
that appeared unchanged for each applied force as the 
overall DNA length was varied from 2.2 kbp to 4.2 kbp, 
whereas the simulation showed a significant increase in 
Az at the longer DNA length. Second, no discontinuity 
was observed in the (directly measured) filtered torque 
data (Fig. 1), yet the simulation predicted a small jump. 

Simulation is not understanding. Here we analyze the 
system theoretically, focusing on the physical causes of 
the behavior at the transition. We use as our framework 
Marko's two-phase coexistence model [6, 10], which we 
generalize to incorporate extra terms that represent the 
interfacial energy between the plectoneme and straight 
regions of the DNA. We show that any model of the 
supercoiling transition in this parameter regime can be 



summarized by four force-dependent parameters. Af- 
ter extracting these parameters directly from the exper- 
iments, we use them to predict the torque jump (which 
we then measure) and to explain why the extension jump 
appears length independent. Finally, we use our formu- 
lation to test various models of plectonemes, finding dis- 
crepancies mainly at small applied force. 

The transition occurs at the critical linking number K* 
when the two states have the same free energy J^, where 

is defined by the ensemble with constant applied force 
and linking number. We therefore need models for the 
free energy and extension z of the SS and CS. 

The properties of stretched, unsupercoiled DNA are 
well-established. At small enough forces and torques 
that avoid both melting and supercoiling, DNA acts as 
a torsional spring with twist elastic constant C [6] [18]: 
J^ss{K, L) = I (27Tf )^ L - FeffL, where K is the added 
linking number, L is the overall (basepair) length of 
DNA, the effective force F^g = F - kT^F/B [6] (see 
supplemental material), F is the force applied to the 
ends of the DNA, B = 43 ± 3 nmx/cT is the DNA's 
bending elastic constant, C — 89 ± 3 nmxA:T, and the 
thermal energy kT = 4.09 pN nm for this experiment 
(at 23.5°C). Differentiating with respect to K gives the 
torque: tss = = ^""^X"- '^^^ extension of un- 

supercoiled DNA is shortened by thermal fluctuations, 
and in the relevant force regime is approximately given 
by ^ss = i{Tss)L, where [11] 
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Since supercoiling theories must include contact forces, 
they are less amenable to traditional theoretical methods. 
Even so, many theories have been successful in predict- 
ing properties of the CS; such methods have included de- 
tailed Monte Carlo simulations [12], descriptions of the 
plectoneme as a simple helix [7, 8, 13], and a more phc- 
nomenological approach [6]. However, none of these the- 
ories has yet been used to predict discontinuities at the 
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FIG. 1: Extension and torque as a function of linking num- 
ber K, for L = 2.2 kbp at _F = 2 pN. Black lines show data 
from Ref. [4], smoothed using a "boxcar" average of nearby 
points. The green lines show worm-like chain (WLC) predic- 
tions below the transition [in the unsupercoiled "stretched" 
state (SS)], and fits to the data after the transition [in the "co- 
existing" state (CS)], linear for the extension and constant for 
the torque. The size of the torque jump, not visible in the 
smoothed data, is implied by the coexisting torque r, the CS 
fit, and the transition linking number K* in the extension 
data. Inset: Simulated DNA showing the CS of a plectoneme 
and straight DNA, ignoring thermal fluctuations. The ends 
are held with fixed orientation and pulled with a constant 
force F, here 2 pN. 

SS-CS transition. Here we connect the free energy and 
extension predictions from any given model to the corre- 
sponding predictions for discontinuities at the transition. 

We will use the framework of two-phase coexistence 
adopted by Marko [6, 10] to describe the CS as consist- 
ing of two phases, each with constant free energy and 
extension per unit length of DNA [19]. Since phase coex- 
istence leads to a linear dependence on K of the fraction 
of plectonemic DNA (keeping the torque fixed), in this 
model both Tcs and zcs £^re linear functions of added 
linking number K and length L (just as the free energy 
of an ice-water mixture is linear in the total energy, and 
the temperature remains fixed, as the ice melts). This lin- 
earity, along with the known properties of the SS, allows 
us to write J^cs and zcs as (see supplemental material) 

Tcs{K,L)=To + 2iTrK-(^^+F,^L- (2) 
zcs{K, L) = -zo - qK + (e(r) -I- ^q) L, (3) 

where q is the slope of extension versus linking number 
and T is the CS torque. That is, Tcs and zcs are specified 
by four force-dependent values: their slopes with respect 
to K (r and q), which describe how the plectonemic phase 
coexists with the stretched phase; and K = L = offsets 



FIG. 2: Directly measuring the torque jump by observing 
thermal hopping, for the same conditions as Fig. 1. As linking 
number K is slowly increased near K* , thermal fluctuations 
induce hopping between states with (CS) and without (SS) 
a plectoneme. Averaging over these two states gives a direct 
way of measuring the torque jump: analogously to a lock- 
in amplifier, we set a threshold in the extension signal to 
separately average the SS (black) and CS (red) data near the 
transition. Using multiple traces, we find an average torque 
jump of Ar = 2.9 ± 0.7 pN nm for L = 2.2 kbp at F = 2 
pN. Additionally, this value of Ar implies (see text) that the 
transition should happen over a range of linking number K 
(top) of about 0.9 turns, as it does. 

{J^o and zq), which describe the extra free energy and 
extension necessary to form the interface between the 
phases — the end loop and tails of the plectoneme. 

The experimental observables can then be written in 
terms of these four values. Easiest are t and q, which are 
directly measured. Next, the linking number K* at the 
transition is found by equating the CS free energy with 
that of the SS: Tcs{K\L) = Tss{K* , L) implies 

r for' 
K*^—{T + AT),withAT=^—To. (4) 

where Ar is the jump in the torque at the transition. 
Lastly, inserting K* from Eq. (4) into Eq. (S12), we find 
the change in extension at the transition: 

Az^zo + q\l^ l(s,{t) - e(r + V^^CT^)) . 

(5) 

To additionally include entropic effects, we can write 
J-Q = ^ — TS, where is the energy cost for the end-loop 
and tails, and S is the entropy coming from fluctuations 
in the location, length, and linking number of the plec- 
toneme. Using an initial calculation of S that includes 
these effects (in preparation; see supplemental material), 
we find that S varies logarithmically with L, and that 
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FIG. 3: The four parameters describing the CS (coexisting 
torque r, extension versus hnking number slope q, and the 
extra free energy To and extension zo necessary to form the 
end loop and tails of the plectoneme), as a function of applied 
force. The circles show values calculated from experimental 
data taken at two different overall DNA lengths L. Model 
predictions for our simulation [20] and Marko's model [6] are 
shown as solid and dashed lines, respectively (using S = for 
To predictions). The circular end-loop model uses average r 
and q values from the experiment to predict To and zo, shown 
as dotted lines. 



setting 5* = is a good approximation except when L 
changes by large factors. 

Given experimental data (r, g, K* , and A2;), we can 
solve for the four CS parameters. The results from 
Ref. [4] are shown as circles in Fig. 3 for the two overall 
DNA lengths tested. If we assume that the DNA is ho- 
mogeneous, we expect the results to be independent of L 
(except for a logarithmic entropic correction to J-q that 
would reduce it at the longer L by about fcTlog2 « 5 
pN nm; see supplemental material). We do expect To 
and zq to be sensitive to the local properties of the DNA 
in the end- loop of the plectoneme, so we suspect that 
the difference in zq between the two measured lengths 
could be due to sequence dependence. With this data, 
we can also predict the length-dependence of the discon- 
tinuities, as shown in Fig. 4 (left). Here we included 
entropic corrections to Tq (see supplemental material), 
and we find that entropic effects significantly decrease 
the length-dependence of the extension jump. 

Note that here we are solving for the experimental 
size of the torque jump using the observed K* and r in 
Eq. (4). We also find direct evidence of At in the data 
by averaging over the torque separately in the SS and CS 
near the transition (Fig. 2) . With data taken at F = 2 pN 
and L — 2.2 kbp, we find At = 2.9 ± 0.7 pN nm, in good 
agreement with the prediction from K* (3.9±2.6 pN nm; 
see Fig. 4). We can also predict the width of the range 
of linking numbers around K* in which hopping between 
the two states is likely (where \/S.T\ < kT): expanding to 
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FIG. 4: (Left) Predicted length-dependence of the extension 
and torque jumps at _F = 2 pN. Using the CS parameters ex- 
tracted from the experiment at two different lengths, Eqs. (4) 
and (5) predict the L-dependence of and Ar. The circles 
show experimentally-measured values [with the torque jump 
here calculated from K* using Eq. (4)]. Without entropic 
corrections to To {S = 0; dot-dashed lines) Az depends no- 
ticeably on L, but including an initial estimate of S (solid 
lines) shows that entropic effects can significantly reduce this 
length-dependence. (Right) Force-dependence of the exten- 
sion and torque jumps, and predictions from two models. 
Disagreements with experimental data can be understood in 
terms of the four CS parameters in Fig. 3. Also plotted as a 
diamond is Ar measured using the direct method depicted in 
Fig. 2. 



first order in K — K* gives a width of 2fcr/(7rAT). This 
predicts a transition region width of about 0.9 turns for 
the conditions in Fig. 2, agreeing well with the data. 

We can now use various plectoneme models to calculate 
the four CS parameters, which in turn give predictions 
for the experimental observables. The results are shown 
as lines in Fig. 3 and Fig. 4 (right). As we expect entropic 
corrections to be small (changing Tq by at most about 5 
pN nm), we set S" = for these comparisons. 

First, we test Marko's phase coexistence model [6]. 
The plectoneme is modeled as a phase with zero exten- 
sion and an effective twist stiffness P < C. Shown as 
dashed lines in Fig. 3, the Marko model predicts the co- 
existing torque and extension slope well, with P as the 
only fit parameter (we use P — 26 nm). However, the 
Marko model (and any model that includes only terms in 
the free energy proportional to L) produces Tq — and 
zo = 0. 

In order to have a discontinuous transition, we must 
include the effects of the end loop and tails of the plec- 
toneme. The simplest model assumes that the coexis- 
tence of stretched and plectonemic DNA requires one 
additional circular loop of DNA. Minimizing the total 
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free energy for this circular end-loop model gives 

JTq = 2n^/2BF^s - 27rrWrioop; (6) 
zo = 2<(t) VB/{2F,s) - gWrioop, (7) 

where Wrioop is the writhe taken up by the loop. For a 

perfect circle, Wrioop = 1, and Wrjoop < 1 for a loop with 
two ends not at the same location. We chose Wrioop = 0.8 
as a reasonable best fit to the data. Using the experimen- 
tally measured t and q, the predictions are shown as solid 
lines in Fig. 3 and Fig. 4; J^q is fit fairly well, but Zq is 
underestimated, especially at small applied forces. 

In an attempt to more accurately model the shape 
of the plectoneme, we use an explicit simulation of an 
elastic rod, with elastic constants set to the known val- 
ues for DNA. We must also include repulsion between 
nearby segments to keep the rod from passing through 
itself. Physically, this repulsion has two causes: screened 
Coulomb interaction of the charged strands and the loss 
of entropy due to limited fluctuations in the plectoneme. 
We use the repulsion free energy derived for the helical 
part of a plectoneme in Ref. [7], modified to a pairwise 
potential form (see supplemental material) . We find that 
the simulation does form plectonemes (inset of Fig. 1), 
and we can extract the four CS parameters, shown as 
solid lines in Fig. 3 [20]. Since JFq and zq are nonzero, 
we find discontinuities in the extension and torque at the 
transition; their magnitudes are plotted in Fig. 4. 

Both the circular loop model and the simulation pro- 
duce torque and extension jumps of the correct magni- 
tude, but in both cases Az has an incorrect dependence 
on force and too much dependence on length. Our ap- 
proach provides intuition about the causes of the dis- 
crepancies by singling out the four values (connected to 
different physical effects) that combine to produce the 
observed behavior. Specifically, we can better under- 
stand why the models' predictions are length-dependent: 
as displayed in Fig. 4 (top left), the negligible length- 
dependence observed in experiment is caused by a sub- 
tle cancellation of a positive length-dependence [smaller 
than either model, and described by Eq. (5)] combined 
with a negative contribution coming from entropic ef- 
fects. One would expect, then, that any plectoneme 
model (even one that explicitly includes entropic fluc- 
tuations) might easily miss this cancellation. In general, 
without this intuition, it is difficult to know where to 
start in improving the DNA models. 

The largest uniform discrepancy happens at small ap- 
plied forces, where both models underestimate zq [21], 
leading to an underestimate of Az. We have exam- 
ined various effects that could alter zo, but none have 
caused better agreement (see also supplemental mate- 
rial). Adding to the circular end- loop model softening or 
kinking [14] at the plectoneme tip, or entropic terms from 
DNA cyclization theories [15, 16], uniformly decreases zq- 
Increasing B in Eq. (6) by a factor of four (perhaps due 



to sequence dependence) does raise zq into the correct 
range, but it also raises !Fo from Eq. (7) to values well 
outside the experimental ranges. Finally, zq would be 
increased if multiple plectonemes form at the transition, 
but we find that the measured values of JFg are too large 
to allow for more than one plectoneme in this experiment. 

Support is acknowledged from NSF Grants DMR- 
0705167 and MCB-0820293, NIH Grant GM059849, and 
the Cornell Nanobiotechnology Center. 
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The language of phase coexistence is approximate in that 
the finite barrier to nucleation in one-dimensional sys- 
tems precludes a true (sharp) phase transition. 
We have also explored increasing the entropic repulsion 
by a constant factor of up to 3. Though this docs bring 
the torques closer to the experiment, the only other sig- 
nificant change is a decrease in J^o (data not shown) — 
specifically, this does not change the discussed discrep- 
ancies between the simulation and experiment. 
[21] Though the 4.2 kbp zo data alone would be arguably 
consistent with the model predictions, the 2.2 kbp data 
highlights the discrepancy at small applied forces. 
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Behavior of extended DNA with fluctuations 

The behavior of extended DNA is appreciably affected by thermal fluctuations. For the applied forces in the range 
considered in this experiment, we can use the following fixed-torque free energy: 
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where the last term is the lowest-order correction due to fiuctuations [11]. 
The fiuctuations decrease the extension: 
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(The —1/32 in Eq. (1) comes from an approximation to a higher-order correction [11].) 
Expanding the last term of Eq. (SI) to match the form of a "zero-temperature" chain, we can instead write 
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where the effective force and twist elastic constant are given by 



F^s = F- kT^ 



Ceff = Chare { I + kT 
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Note that Cefi is a function of force: there is less "softening" at higher forces. In the experiments of Forth et al., 
the renormalized Cgs was measured directly via the torque. However, the range of applied forces was small enough 
that Coil did not change appreciably, and a single value of C = (89 nm)A;T was quoted. Here, we also use the same 
renormalized but force-independent value for C. 

Changing Eq. (S3) to a fixed-linking-number expression via a Legendre transformation, we arrive at our expression 
for the straight state free energy (also found in Ref. [6]): 



C / K 
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Derivation of linear expressions for J^cs and zcs 

We first write down the linear scaling of the free energy and extension with linking number. For any SK that does 
not take the system out of the CS, 

J^cs{K + SK, L) = Tcs{K, L) + 2t,t5K- (S7) 
^cs(^ + SK, L) = zcs{K, L) - q6K, (S8) 

where q is the slope of extension versus linking number and r is the CS torque. Next, to find the scaling with 

increasing L, we imagine adding a piece of stretched DNA of length SL at the coexisting torque (keeping the system 
in a stable CS). This also adds an amount of linking number that scales with 5L, dK[SL] = t5L/{2ixC), which we 
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will have to unwind to get back to the original K. First adding the piece of stretched DNA, and then unwinding to 
find the dependence on L only, we find 

J'csiK, L + SL)= Tcs{K, L) - + F,^ 5L. (S9) 

Similarly for the extension, [using ^(r) from Eq. (1)] 

zcs{K,L + SL) = zcs{K,L) + (^(r) + ^,q) 5L. (SIO) 

Combining Eqs. (S7) and (S8) with Eqs. (S9) and (SIO), we can write the free energy and extension of the CS as 
linear in K and L, each with a slope and an intercept: 

J'csiK, L)=To + 27TTK - + i^eff) L- (Sll) 

zcs{K, L) = -zo -qK+ (^(r) + L. (S12) 

Note that C and ^(r) are known from experiments on stretched DNA, leaving the four anticipated force-dependent 
quantities to be described by a theory of supercoiling: r, q, Tq, and zq. 



Self-repulsion 

It is essential to include a repulsive force between sections of the DNA that come near each other; without it, the 
rod can pass through itself, unphysically removing linking number in the process and preventing the formation of 
plectonemes. The physical origins of repulsive forces in DNA include both electrostatic and entropic effects. We use 
discretized versions of the repulsive interactions described in Ref . [7] . 

Electrostatic forces are modeled using a Debye-Huckel screened Coulomb interaction: 

Esc{r) = J ^ , S13 

€ r 

where v = 8.4 nm~^ is the effective number of electron charges per unit length, Xd = 0.8 nm is the Debye screening 
length, and e?./e = 2.9 pN nm^. (These values are dependent on the ionic concentration of the buffer, and were picked 
to match with « 150 mM NaCl.) 

The entropic free energy of a helical structure is calculated in Ref. [7] , coming from the increasing confinement of 
fluctuations in more tightly coiled structures. We use the same free energy, written as a pairwise interaction between 
segments: 

_ 25/3V^r(l/3) kTd^ 

Since we also include straight parts of the DNA that should not have the same entropic interaction, we cut off the 
entropic potential at a distance of 2B/kT, where the argument for the form of the potential breaks down [7]. 



Extra terms in the circular end-loop model 

Extra terms in the free energy that we have not considered would change the predictions of the circular end-loop 
model — these could include electrostatic interactions, entropic effects, etc. In fact, we can solve for the properties that 
such an extra free energy term (call it .Fextra) would need to have in order to make the model match the experimental 
data. 

Adding this unknown term, we have 

J^i{Ki, Li) = ^ [2t,{Ki - Wrioop)]' + (STr)^ A + T,^,,^{Ki, Li). (S15) 
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FIG. SI: Entropic corrections from the literature do not help the circular end-loop model fit the data. The dots show the 
required free energy contribution .7-"extra (top plot) and its derivative with respect to end-loop circumference dTcyLtT&l dLi (bottom 
plot) that would produce an To and zo that match with the experiment (with Wrioop = 0.8). Bars on the top plot show the 
required derivative, the value of which is shown on the bottom plot. Vertical grey lines show one standard deviation error 
bars. Note especially the inability of any of the proposed entropic terms to match the well-constrained negative derivative at 
large end-loop circumferences (which happen at low force in the experiment); this produces Li (and thus zo) that are too small 
at low forces. A lessening of the effective force felt by the end-loop of about 0.5 pN would help agreement, but none of the 
proposed corrections provides this. 

Since the terms we will imagine adding will not depend on Ki, we will assume that Tcxtxa is only a function of Li. 
Setting the force and torque equal to the coexisting state values {dJ-i/dLi ~ — (i^cff + t^/(2C)); dTi/dKi = 27rT) then 
gives 
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We now use the fact that 
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to solve for the necessary values of ^extra and dJ^e:x.tia/dLi in order to match with the experimental J^o and zq. We 
find 



^extra = J^O + 27rrWrioop - i^effi;* - -jr^ (S20) 



dLi L* 

where 



- i^eff, (S21) 

These required properties of the added free energy term are plotted in Fig. SI for Wrioop = 0.8. 

Wc can then test whether different possible extra free energy terms would match the requirements. liere we try 
four possibilities taken from the literature. First, there is electrostatic repulsion coming from like charges on opposite 
sides of the DNA circle. This looks like (using the Debye-Huckel formulation from Ref. [7]) 

_^^ectrostatic ^ kTlg^^Ko ( ^] h (S23) 



ttAd 

and is plotted in yellow in Fig. SI. Second, Odijk calculates the free energy for a circular DNA loop and finds terms 
in the free energy [16] [Eq. (2.13)] 

« = ^ri„g^-m!.; ,S24) 

this is plotted in purple in Fig. SI. Third, a similar term is found by Tkachenko in solving for the J-factor for 
unconstrained DNA cyclization [18] [Eq. (4)]: 

^ex'^?f = 5fcT log ; (S25) 

this is plotted in green in Fig. SI. Finally, we could imagine that entropic contributions from confinement similar 
to the one used by us for our elastic simulation could be important. Although the form was derived for a different 
configuration (superhelical DNA) , we could try it to see if something similar might help. Integrating the confinement 
entropy from Marko and Siggia [7] over a circle gives 

kT 

^confinement r ^^9^'! 

-^extra - (B/fcT)l/3(i/(27r))2/3 ^' 

which is plotted in blue in Fig. SI. 

Although these possible terms are only initial guesses at the possible corrections due to entropic and other effects, 
we see that they are all qualitatively unable to help, especially at long loop lengths, which is where the circular loop 
model fares worst at fitting the data. 



Calculating entropic contributions from fluctuations in plectoneme location, length, and linking number 

To investigate entropic effects, we would like to find the free energy of states with multiple plectonemes [19], 
including fluctuations of linking number and length both within individual plectonemes and moving among different 
plectonemes. We can achieve this by calculating the partition function for a state with n plectonemes, identifying 
unique states by the plectoneme positions Sj, the plectoneme lengths Lpj, and the plectoneme linking numbers Kpi: 



Zn{K,L) = 




(S27) 



exp [-J^ri{L, K, Lp^, Kpi)/kT], 
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where we have neglected the compHcations coming from the possibihty that plectonemes could overlap. The constants 
Lq and Kq set the length change and linking number change, respectively, that produce an independent state. Since 
we are only concerned with the free energy difference between the straight state and coexisting state, these constants 
would be set by the change in entropy of the degrees of freedom in the straight state that axe lost to the collective 
modes we are integrating over in the coexisting state. 

The first line of integrals represents the choice of where to put each plectoneme, which does not change the free 
energy does not depend on Si). We therefore simply get a factor of L", divided by n\ since plectonemes are 
indistinguishable: 



Zr,{K, L) = ^^dMl /'"' TT dLpi r n dKpi exp [ - ^„(L, K, {L^,}, {Kpi))/kT] . (S28) 

n\ J^qJ^o Jo ~ J-oo 

Next we need to know the free energy of coexisting states that are away from the equilibrium plectoneme length 
and linking number. Assuming that the plectoneme free energy density is quadratic in linking number density (as in 
Marko's model [6]), this turns out to be 



(L - ^ Lpi) - Feff{L - ^ Lpi) + nil, 



(S29) 
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where /x is the chemical potential for plectoneme ends and v = 2CFeff/T^. 

We first evaluate the integrals over Kpi, which amount to n Gaussian integrals; this gives 



WdLpi 



Hi Lpi/ci 



1/2 



1 + {1 + v)tS±^^ 
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exp(- 



kT 



{2TrKy 



L-Y.Lpi + {l + v){Y.Lpi) 



-r - Feff (L - ^ Lpi) + niJ. 



)• 



Now changing to unitless variables Xi = Lpi/Lp and y = Lp/L, and rearranging to move all the factors that depend 
on the sum of the plectoneme lengths y into the exponent, the term in the exponent becomes 



kT\ l + vy 



Feffi(l - y) + n/i 1 + - log 



l + vy 



(S31) 



and we have 
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n! L'^K^ 



^ dLp 6 (J2 Lp^ -Lp) dLpi^jllLp^/c^cxp [-f{Lp/L)] 

exp [-f{Lp/L)] 



TT ' 



^0 



' rfy exp[-(/(y) — logy)]. 



The integral in large square brackets (characterizing fiuctuations in the individual plectoneme lengths that do not 
change the total plectoneme length) gives a numerical constant 7„ = 7r"/^/(2"r(3n/2)) = 2^^^W^^^ /{2in — 2)!!. To 
evaluate the y integral over total plectoneme length, we make a Gaussian approximation [noting that the total length 
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is well-constrained by f{y)]. Then the fluctuations in the (fractional) total length of plectonemic DNA are of size 

-1/2 



3n-2 



logy 



(S33) 



where y* is the equilibrium value of y, and the derivative is 



3n- 2 



logy 



2 (1 - y)2 
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LkT{l+vy) I 



y^ 



(1 + vyy 
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Without the entropic corrections, the equilibrium length is y* = {u — l)/v, where u = 2'!tCK / [tL). We can safely 
use this value if we are far from y* =Q and y* = 1, and get 



<y-u = 



( 1 2t'^L 1 



- ^ + 



+ 



3n- 2 



V \ukTC V? {v-u+iy (u-l)2 



-1/2 



(S35) 



[Since we are usually near y* = at the transition, to calculate the length-dependence shown in Fig. 4 (left), we 
approximate y* numerically and use Eq. (S33) instead of Eq. (S35) .] In the end, we have 



u-l 



(3ra-2)/2 



V — U+1 



UV 



1/2 



The full partition function for all plectonemic states is then 

oo 

Z{K,L) = Y,Zr,{K,L) 



eyi^[-T{K,L)/kT]. (S36) 



(S37) 



(which we can numerically approximate by truncating the series at a reasonable n), such that the coexisting state free 
energy is given by Tcs{K, L) = — fcTlog Z{K, L). For the experimental values, we find that only the single plectoneme 
n = 1 state contributes significantly near the transition. 



Independence of results on entropic effects 

In the paper, we have set the entropy from the previous section to zero (5 = 0) for most of the calculations. How 

would wc expect that including S would change any of the results? 

First, S would create a shift between the experimental J-q and the predictions from models that do not include 
fluctuations. We find that this shift is largely independent of force, and is mostly dependent on Lq. We do not 
currently have a way of calculating Lq , but wc expect that it should be on the order of the persistence length of DNA, 
about 50 nm. We find that setting Lq to about 100 nm makes the prefactor equal to 1, or equivalently sets S = Q. If 
we assume that Lq is about equal to the persistence length of DNA, we expect that we would need to shift the model 
predictions by at most about fcTlog2 w 5 pN nm. 

Second, we find that S has a logarithmic dependence on L. This means that we expect JFq to decrease by something 
on the order of kT\og{L2/Li) when we increase the length from Li to L2. For the experimental lengths (with 
L2 ~ 2Li), this again corresponds to a shift of about 5 pN nm. 

Shifting J^Q by these amounts would slightly change only the theory curves for J^q (about 5 pN nm), (about 10 
nm), and At (about 1 pN nm). 
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